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The thermal history of the universe before the epoch of nucleosynthesis is unknown. The maxi- 
mum temperature in the radiation-dominated era, which we will refer to as the reheat temperature, 
may have been as low as 0.7 MeV. In this paper we show that a low reheat temperature has impor- 
tant implications for many topics in cosmology. We show that weakly interacting massive particles 
(WIMPs) may be produced even if the reheat temperature is much smaller than the freeze-out 
temperature of the WIMP, and that the dependence of the present abundance on the mass and the 
annihilation cross section of the WIMP differs drastically from familiar results. We revisit predic- 
tions of the relic abundance and resulting model constraints of supersymmetric dark matter, axions, 
massive neutrinos, and other dark matter candidates, nucleosynthesis constraints on decaying par- 
ticles, and leptogenesis by decay of superheavy particles. We find that the allowed parameter space 
of supersymmetric models is altered, removing the usual bounds on the mass spectrum; the cos- 
mological bound on massive neutrinos is drastically changed, ruling out Dirac (Majorana) neutrino 
masses m„ only in the range 33 keV ~ m„ < 6 (5) MeV, which is significantly smaller from the the 
standard disallowed range 94 eV IS m„ < 2 GeV (this implies that massive neutrinos may still play 
the role of either warm or cold dark matter); the cosmological upper bound on the Peccei-Quinn 
scale may be significantly increased to 10 16 GeV from the usually cited limit of about 10 12 GeV; and 
that efficient out-of-equilibrium GUT baryogenesis and/or leptogenesis can take place even if the 
reheat temperature is much smaller than the mass of the decaying superheavy particle. 
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I. INTRODUCTION 



The initiation of the radiation-dominated era of the universe is believed to result from the decay of coherent 
oscillations of a scalar field whose energy dominated the universe before decay. The decay of the coherent oscillations 
of the scalar field and the subsequent thermalization of the decay products is known as reheating.^ 

The reheat process is often associated with the final stage of inflation. However, reheating could have been episodic, 
with several reheat events after inflation. We will be interested in the final reheating before primordial nucleosynthesis, 
which may just as well have been the result of the decay of a weakly coupled scalar field unrelated to inflation, for 
instance a modulus. For this reason the scalar field </>, whose decay leads to reheating, will not be referred to as the 
inflaton. 

A common assumption is that many of the interesting cosmological phenomena accessible to present-day obser- 
vations occurred after reheating, during the radiation-dominated phase of the early universe. This seems to be a 
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1 While we discuss reheating as the decay of coherent field oscillations, we only make use of the fact that the energy density of 
coherent field oscillations scales in expansion as a~ 3 , where a is the Friedmann-Robertson- Walker scale factor. One could just 
as easily imagine that the universe is dominated by some unstable massive particle species, rather than coherent oscillations of 
a scalar field. 



reasonable assumption, since inflation [jjj erases any initial condition on the number densities of ordinary particles, 
while the reheat process repopulates the universe. The assumption of an initial condition of thermal and chemical 
equilibrium in a radiation-dominated universe is then equivalent to the hypothesis that the maximum temperature 
obtained during the radiation-dominated era, Trh , is larger than characteristic temperatures of cosmological processes 
under investigation. 

The fact that we have no physical evidence of the radiation-dominated era before the epoch of nucleosynthesis (i.e., 
temperatures above about 1 MeV) is a simple, but crucial, point. Therefore, a priori one should consider Trh as an 
unknown quantity that can take any value as low as about 1 MeV. Indeed, there are good physical motivations for 
studying cosmologies with very low Trh ■ 

In theories in which the weakness of gravity is explained by large compactified extra dimensions || , emission of 
Kaluza-Klein gravitons in the bulk constrains the normalcy temperature at which the radius of the compactified 
dimensions is stabilized with vanishing energy density in the compactified space to be in the MeV to GeV region 
||. In practice, the stabilization occurs in a post-inflationary phase, and therefore we can identify the normalcy 
temperature with Trh ■ 

Supergravity and superstring theories usually have particles, such as a gravitino or a modulus, with only gravita- 
tional interactions. Late decay of these particles may jeopardize the success of standard big-bangnucleosynthesis [Q. 
This problem can be solved by assuming sufficiently low Trh, of the order of 10 8 to 10 10 GeV ||. Moreover, it has 
been recently realized || that nonthermal production of these gravitational relics during the inflationary phase can 
impose upper bounds on Trh as low as 100 GeV. 

In this paper, we will show that the phenomenological point of view that the reheat temperature may be as low as 
1 MeV has rich implications for particle dark matter, neutrino mass limits, axion cosmology, and baryogenesis. 

The key point of our considerations is that reheating is not an instantaneous process. On the contrary, the radiation- 
dominated phase follows a prolonged stage of matter domination during which the energy density of the universe is 
dominated by the coherent oscillations of the field <f>. The oscillations start at time HJ 1 and end when the age of the 
universe becomes of order of the lifetime ri 1 of the scalar field. At times HJ 1 < t < ri 1 the dynamics of the system 
is quite involved. During this stage the energy density per comoving volume of the <fi field decreases as exp(— T^t) 
and the light decay products of the scalar field thermalize. The temperature T of this hot plasma, however, does not 
scale asTaa -1 as in the ordinary radiation-dominated phase (a is the Friedmann-Robertson- Walker scale factor) 

, but reaches a maximum Tmax ~ (HrMpi^T^H (Mpi is the Planck mass) and then decreases as T oc a 3 ' , 
signalling the continuous release of entropy from the decays of the scalar field. This scaling continues until the time 
t ~ ri 1 when the radiation-dominated phase commences with temperature Trh- Therefore, before reheating is 
completed, for a given temperature the universe expands faster than in the radiation-dominated phase. Notice that 
Trh is not the maximum temperature during the reheat process. On the contrary, Tmax can be much larger than 
Trh- The behavior of the universe during reheating is discussed in detail in Sec. [n| 

In Sec. Ill we will use these results to compute the relic abundance of a dark-matter species (X) produced during 
reheating. We will consider the case that Trh is smaller than the X freeze-out temperature]^] Although naively 
one might expect a negligible X number density under these circumstances, we find that the X relic density can 
reach cosmologically interesting values. We also show that because of entropy release during the reheat stage, for a 
given mass and cross section the present X abundance is smaller than obtained assuming freeze out in the radiation- 
dominated era. This relaxes the bounds coming from requiring that flxh 2 < 1. Moreover, the parametric dependence 
on the mass and the annihilation cross section of the present abundance is nonstandard. Therefore, most of the 
cosmological constraints on specific particle properties have to be revisited. 

The value of the maximum temperature during reheating, Tmax, delineates different regions for our results. If 
Tmax is smaller than the X mass, X particles generated by collisions in the thermal bath during the coherent (\> 
oscillations are always nonrelativistic. If the cross section is small enough, the X particles do not reach chemical 
equilibrium and the present abundance is proportional to the cross section, Qx oc (a a \v\), in contrast with the usual 
radiation-dominated case in which fix oc (<7a|^ Avoiding overclosure of the universe imposes a lower bound on 
the annihilation cross section. On the other hand, if the X particles reach chemical equilibrium before reheating is 
completed, the ratio of the X number density to entropy density in a comoving volume does not stop decreasing when 
the particles freeze out because entropy is released until reheating is over. Therefore, the present abundance does 
not depend only on (ctaM) -1 (as in the standard case), but also on the reheat temperature Trh- the lower Trh, 
the smaller the abundance. We will provide a formula for fix which reproduces both the standard result when Trh 
is equal to the freeze-out temperature of the relic particle and the nonstandard result Qx oc (caM) when the value 



2 For a different perspective, see Ref. fb 



cross section becomes smaller than some critical value. 

If Tmax and Trh are larger than the X mass, we will show that the relevant processes determining the present X 
abundance occur during reheating or afterwards. 

An important result is that Tmax (or equivalently Hi) is relevant when deciding if X is relativistic or not, but does 
not appear in the final expression for the relic abundance. Therefore, once Tmax has determined the pertinent case, 
fix depends on the physics of the <f> field only through Trh- It is easy to understand why. The X number density 
results from the competition of two rates, the interaction rate and the expansion rate of the universe. Before reheating 
is completed the expansion rate depends only on Tpa , H ~ T /T^ H Mpi, and therefore the final abundance depends 
only upon T RH . 



In Sec. IV we discuss the applications of our findings to some popular cold dark matter candidates. Here, we 
preview some of our results. 

While excluded in the usual scenario, thermal WIMPs with mass larger than the unitarity bound of a few hundred 
TeV Jll|] may be viable dark matter candidates in low reheat models. 

Of the many WIMP candidates, the best motivated seems to be the neutralino, the lightest supersymmetric particle 
(LSP) |Q. In the case in which the LSP is mainly a Bino, requiring Q,b ~ 1 gives an upper bound on the slepton 
mass fri£ R ]l3|-p"5||. However, once we relax the assumption that the reheat temperature is higher than the freeze-out 
temperature of the WIMP, we will show that the upper bound on rhe R is drastically relaxed and may completely 
disappear. The same argument can be applied to other supersymmetric candidates. 

Another striking application for dark matter that illustrates our point is the computation of the relic abundance of 
massive neutrinos. The well-known cosmological Cowsik-McClelland — Lee- Weinberg bound [ fujlrijl rules out neutrinos 
more massive than roughly 90 eV and lighter than around 2 GeV. This result has a significant impact in cosmology, 
ruling out, for instance, the possibility that neutrinos are warm dark matter. This standard result, however, assumes 
that the reheat temperature is much higher than 1 MeV and that neutrinos have been in chemical equilibrium up to 
temperatures of the order of 1 MeV. We will show that if the reheat temperature is as small as allowed by big-bang 
nucleosynthesis, then massive, stable Dirac (Majorana) neutrinos are compatible with cosmology if they are lighter 
than about 33 keV or heavier than about 6 (5) MeV. This implies, for instance, that neutrinos may still be warm or 
cold dark matter and play a significant role in other cosmological or astrophysical phenomena. 

We then proceed by investigating the implications of a low reheat temperature for axion cosmology. It is well 
known that in the standard scenario the oscillations of the axion field generated by the misalignment mechanism 
overdose the universe unless the Peccei-Quinn scale, fpQ, is smaller than about 10 12 GeV. This bound, however, is 
obtained assuming that the reheat temperature is larger than the QCD scale. When this assumption is abandoned, 
the cosmological upper bound on fpQ is significantly relaxed to Jpq £ 10 16 GeV if Trh ~ 1 MeV. 

If the X particle has a nonvanishing but small decay lifetime tx , its decay products may destroy the light elements 
generated during primordial nucleosynthesis. This gives strong constraints in the plane (Mx,tx) |l§| ] that are very 
sensitive to the number density of the species X at freeze out. Lowering the reheat temperature implies a smaller 
number density and therefore much less restrictive bounds. 

These results all imply that presently stated cosmological limits may not always be relevant in limiting particle 
properties such as the supersymmetric mass spectra in the experimentally verifiable range of future colliders. 

As a last application, in Sec. ^ we analyze the production of unstable superheavy states during the process of 
reheating, keeping in mind the possibility that the subsequent decay of these states may generate the observed baryon 
asymmetry. The fact that Tmax is larger than the reheat temperature may give rise to an efficient production of 
these superheavy states. As a result, out-of-equilibrium GUT baryogenesis and/or leptogenesis can take place even if 
the reheat temperature is much smaller than the mass of the superheavy decaying particle. This is particularly useful 
in supersymmetric scenarios where Trh has to be low enough to avoid the overproduction of gravitinos and other 
dangerous relics. 



Finally, Sec. VI summarizes our results. 



II. THE DYNAMICS OF REHEATING 

A. The Relevant Boltzmann Equations 

In this section we study the Boltzmann equations for the time evolution of a system whose energy density is in the 
form of unstable massive particles 0, stable massive particles X, and radiation R (other similar studies can be found 
in Refs 
and anni 



9[-pT[ ) . We assume that <f> decays into radiation with a rate T^, and that the X particles are created 
rilate into radiation with a thermal-averaged cross section times velocity (crv). The corresponding energy 



and number densities satisfy the differential equations 
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Here, we assume that each X has energy (E x ) ^ VM 2 + 9T 2 and the factor 2(E X ) is the average energy released in 
X annihilation. Later we will assume px — (Ex)nx- The Hubble expansion parameter H is given by 
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8tt 
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The equilibrium number density for particles obeying Maxwell-Boltzmann statistics can be expressed in terms of K2, 
the modified Bessel function of the second kind: 
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where g is the number of degrees of freedom of the X-particle species. 

For cosmological considerations it is more appropriate to express in terms of the reheat temperature T R h using 
the conventional expression 
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where g*(T) describes the effective number of degrees of freedom at temperature T. This expression defines T R h- 

Next, we express Eqs. in terms of dimensionless variables and convert time derivatives to derivatives with 

respect to the scale factor a. The dimensionless variables we choose are 
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RH 



'a 3 



R — prcl 4 ; X = nxa 3 ; A = a/aj 



(7) 



The choice of in the definition of 4> is for convenience; any mass scale would suffice. The factor aj will be chosen 
as the initial value of the scale factor for the integration. Since no physical result can depend upon the choice of aj, 
we are free to choose 
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In terms of the new variables, Eqs. 
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At early times the energy density of the universe is completely dominated by the <f> field. The initial value of the <fi 
energy density can be expressed in terms of the initial expansion rate, f/j, as p^ = (3/8ir)Mp,Hj. Therefore, we will 
solve Eqs. (|9|)~(jll|) choosing the following initial conditions: 
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B. The Temperature— Scale Factor Relation 



During the epoch between the initial time, Hj , and the completion of reheating at time T^ 1 , the temperature of 
the universe does not scale as T ~ a -1 as in the radiation-dominated era, but follows a different law Q|. This unusual 
relation between the temperature and the scale factor, derived below, will significantly affect the calculation of the 
relic abundance of X particles. 

The temperature of the system is measured by the radiation energy density, and therefore T is related to R by 
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At early times (H ^> r^), we can approximate the right-hand side of Eq. ( |To| ) by retaining only the terms proportional 
to the <f> energy density and by taking $ ~ $/. The solution of Eq. ([l0|) then becomes 



Using this result in Eq. we obtain the expression of the temperature T as a function of the scale factor, 

T = T MAX f(A) , 

where Tmax and f(A) are given by 
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The constant k is defined as 
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The function f(A) starts as zero, then grows until A$ — (8/3) 2 / 5 , where it reaches its maximum f(Ao) — 1 (cor- 
responding to T = T M ax), an d then decreases as A -3 / 8 . Therefore, for A > A , Eq. ( |l5| ) can be approximated 
by 
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This result shows that during the phase before reheating the temperature has a less steep dependence on the scale 
factor than in the radiation-dominated era. In other words, as the temperature decreases, the universe expands faster 
before reheating than in the radiation-dominated epoch. Notice also that Tmax can be much larger than Trh, as 
long as Hi > T RH /Mpi. 



C. The Temperature— Expansion Rate Relation 



Next, consider the temperature dependence of the expansion rate H during the epoch of reheating. Between the 
time when Tmax is obtained and the decay time T^" 1 , the scalar field energy density scales as — ^jT RH A~ 3 . Since 
H 2 ~ (8tt/3)p4 > /M P[ , we can express H as 



r2 = 8tt $iT rh A- 
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Now we can use Eq. (03) to express A in terms of T, with the result 
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This result can be compared to the result for a radiation-dominated universe (H oc T 2 ) and a matter-dominated 
universe 



III. CALCULATION OF THE RELIC ABUNDANCE 



In this paper we are interested in the situation in which the X particles never obtain chemical equilibrium in the 
radiation-dominated era after reheating.^ This means that Trh must be smaller than the conventional freeze-out 
temperature (roughly equal to Mx/20 for nonrelativistic, weakly interacting particles). In this situation we can 
encounter several possibilities. In the first case, the X particles are always nonrelativistic and never in chemical 
equilibrium, either before or after reheating. In the second case, the nonrelativistic X particles reach chemical 
equilibrium, but then freeze out before the completion of the reheat process. Finally, we can consider the case when 
the relevant processes of particle production and freeze out occur before reheating at temperatures at which X is still 
relativistic. In this section, we will compute the X thermal relic abundance in all these cases. 

In principle, the relic abundance could receive contributions from other sources, like the direct (f> decay into X 
particles p2| , or from the production and decay of heavy particles eventually decaying into X. In this paper we 
will ignore these model-dependent effects and therefore our calculation can be viewed as a lower bound on the X 
abundance. 

The nonrelativistic and relativistic cases are discriminated by the conditions Mx > Tmax and Mx < Imax, 
respectively. This translates into a condition on if/; for instance, the nonrelativistic case corresponds to 
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In this paper we will treat Trh and Tmax as free parameters and we will not rely on particular models of inflation 
or (f> decay. Nevertheless, it is useful to show what kind of cj> physics can give rise to the different cases considered in 
this section. First of all, we are interested in very low reheat temperatures. This can be achieved if (f> has a typical 
gravitational decay width ~ M?/Mj, ; , for which 
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The value of Tmax is determined by the initial <p energy density. If Ma, is the mass scale characterizing the physics 
of 4>, we can expect p,p(ai) ~ Mf . This happens, for instance, in hybrid models of inflation. In this case we find 
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On the other hand, in the case of chaotic inflation, one expects that the <fi field has an initial value of the order of the 
Planck mass and therefore p • 4,(0,1) ~ M?Mj> ; . This leads to 
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For an X particle with typical electroweak mass and for Ma, ~ 100 TeV, the two options correspond to the nonrela- 
tivistic and relativistic case, respectively. 



3 Here we make the usual distinction between chemical equilibrium and kinetic equilibrium. Kinetic equilibrium can be achieved 

by JVnumber conserving scatterings, such as 7V < > jX (7 represents a light degree of freedom like a photon). Chemical 

equilibrium can only be achieved by processes that change the number of X particles, such as XX < > 77. For massive 

particles the cross section for the second process may be orders of magnitude smaller than the cross section for the first process, 
and it is possible to assume kinetic equilibrium but not chemical equilibrium. 



A. Nonrelativistic Nonequilibrium Production and Freeze Out 



Let us suppose that the X particles are always nonrelativistic and the condition in Eq. ( |2l| ) is satisfied. Since we 
are considering the case in which X does not reach chemical equilibrium (X -C X eq ), at early times, Eq. ( pd| ) can be 
approximated by 

S = (^) 1/2 A ~*** ^WpiThhXI^ . (25) 
The equilibrium distribution in the nonrelativistic limit is 

Xe " = gA * (2^:) ' e M~Mx/T) , (26) 



where g is the number of degrees of freedom of the particle X and the temperature T is given by Eq. ( |18| ) . We express 
the thermal-averaged annihilation cross section times velocity as 

M S lk (°" + ■ (27) 

Here the dimensionless coefficients a s and a p describe, respectively, the s-wave and p-wave annihilations in a nonrel- 
ativistic expansion of the cross section. Using Eqs. fll2]), (|l8|), (p6|), and ( p7| ) in Eq. (p5|), we obtain 
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We will find the solution for X(oo) as a Gaussian-integral approximation to Eq. (|28|). Although Eq. (28) is valid 
only at early times, we will integrate it in the full range between A = and A = 00. This is a good approximation 
because the exponential suppression makes the right-hand side of Eq. negligible anywhere outside a small interval 
of scale factors centered around A = A* , with 

TTm7"J for,- wave, (29) 
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Using Eq. (|T^), we find that A* corresponds to a temperature T» = AMx/17 for s-wave and T„ = 4Mx/l5 for 
p-wave. Therefore T» is the temperature at which most of the A-particle production takes place. We will assume 
T* < ?max or else the final A-particle density is suppressed by a very small exponential function. The Gaussian- 
integral approximation to Eq. (Eq) isn 
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Next we want to relate Xoo to the mass density of X particles today. After particle production stops at A ~ A*, 
the factor X oc nxA 3 = Xoo remains constant. Therefore, at reheating 

Px{T RH ) = M x nx{T RH ) = MxX^T^A^ , (32) 

where from Eq. ( |T8] ) 



4 The result of the Gaussian integration in Eq. ( pi] ) is a very good approximation of the exact integral, which is 
given by = {%g 2 M x ) / (?>-k 3 H 1) [a s f s (2Mx/ K T mal ) + a p f p (2Mx / K,T ma J)], where f s (x) = exp(~x)8\J2l =0 x k ~ 12 /k\ ~ 
v / 2^(17/2) 17/2 exp(-17/2)x" 12 and f p {x) = exp(-£c)(8!/4) ^Lo xk ~ 12 ~ 2^ (15/2) 15/2 exp(-15/2)a:- 12 ~ / s (a;)/4. 
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Also, at reheating the radiation energy density is 



Pr(Trh) = 7^9*(T RH )T RH . 
After the completion of reheating the universe is radiation dominated, and 
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Of course the extraction of energy from the scalar field is not an instantaneous process, but we can use T = Trh in 
Eq. ( |35| ) and correct for the entropy release after Trh. It is straightforward to demonstrate (see the Appendix) that 
relatively independent of the model parameters, only about 25% of the comoving <j> energy density has been extracted 
at T = Trh- At temperatures smaller than Trh some residual entropy is released by the decays of the scalar field till 
the time when the energy density in radiation significantly dominates over the energy density of the scalar field. One 
can show that (again nearly independent of model parameters) the comoving entropy increases by about a factor of 
8 after Trh ■ Therefore all th e ana lytic estimates should be divided by a factor of 8. This is confirmed by numerical 
calculations as shown in Sec. HIE. The reader is referred to the Appendix for more details. 

From Eq. (|3^) with T = Trh and the extra factor of 1/8, we obtain an estimate for the present energy density of 
X particles in units of the critical density 
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(nonrelativistic nonequilibrium production during reheating era) . 

Here we have used T now = 2.35 x 10~ 13 GeV and ^IrIt? = 4.17 x 10~ 5 , including the contributions from the cosmic 
background radiation and from neutrinos. 

Notice that in this case fix is proportional to the annihilation cross section, instead of being inversely proportional, 
as in the case of the usual thermal relic abundance calculation in a radiation-dominated universe. 

The basic assumption used in this section is that the X particles never reach thermal equilibrium. This hypothesis 
holds if at the time of maximum particle production the X number density is less than the equilibrium value, 
Xqo < A eg (T , »). Using Eqs. (|26| ) and (|3l|), we find that this condition corresponds to a limit on the annihilation cross 
section a s < a s for s-wave or a p < a p for p-wave, with 
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In the non-equilibrium case considered in this section, we find that there is a maximum value of Qx that can be 
achieved. This is obtained by replacing the constraint a s < a s into Eq. (|3C" 
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A similar constraint can be obtained in the case of p-wave annihilation. 



To conclude this section, we want to show that for the case under consideration the process of particle production 
always freezes out before reheating, i.e., 

H(A*) > T . (40) 

Since the universe is matter dominated by the <p field, H scales like A -3 / 2 and therefore H(A st ) = HiA* 3 ^ 2 . Using 
Eqs. (p9|)-(^0|), we find that the relation in Eq. (^) is satisfied whenever Mx > "ATrh- This condition is always 
verified once we assume that the X particles do not thermalize after reheating. 



B. Nonrelativistic Equilibrium Production and Freeze Out 



We will now consider the case in which the annihilation cross section is large {a s > a s or Op > a p ) and the X particle 
species reaches chemical equilibrium before reheating (this case is also discussed in Ref. 19 1 ) . The calculation of fix 
is now analogous to the ordinary calculation of thermal relic abundances. However, the result is different because 
the relation between temperature and scale factor is not the same as in the ordinary case of a radiation-dominated 
universe. 

The freeze-out temperature Tp is obtained by solving for the condition 



n^{T F ){av) 



H(T F 



(41) 



Here = X eq A 3 T^ H is the equilibrium number density of X particles and the expansion rate H as a function of 



temperature was expressed in Eq. (Bfj). 



Defining xp = Mx/Tp, the condition in Eq. (41) can be written as 
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(42) 



The factor 5/4 in front of the p-wave term has been added to match the analytic solution of the Boltzmann equation. 
Notice that Eq. ( (42| ) admits a solution only for a s > 1.5a s or a p > 0.4a p , in nice agreement with the starting 
assumption on the annihilation cross section. 

For comparison, we remind that in the case of a radiation-dominated universe the expansion rate is H = 
V / 47r 3 .g !i ,/45r 2 /Mpi and the analog to Eq. (|l|) is 
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In order to compute fix we can use Eq. ( |3q ) (again with the understanding that there will be an overall correction 
due to the fact that the reheating process is not complete at Trh) to derive the result 
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Thus we obtain, in the case decoupling occurs during reheating, the result 
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(nonrelativistic equilibrium production during reheat era) . 

In this case fix is inversely proportional to the annihilation cross section, as in the radiation-dominated case. Eq. 
( [45| ) generalizes the ordinary result of decoupling during the radiation-dominated era, which is given by 
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(nonrelativistic production during radiation era) 



Indeed, Eq. (^|) approximately reduces to Eq. ( phf ) as Trh approaches Tp. Furthermore, Eq. (^h|) also reproduces 
Eq. ( pq ) when a s approaches a s from above. 

The effect of a low reheat temperature is to reduce the relic abundance with respect to the ordinary case by a 
factor T RH Tp d / (T£ ew ) 4 , where T F M and T£ ew are the freeze-out temperatures in cosmologies with high and low Trh , 
respectively. This suppression factor can be understood in the following way. During the epoch before reheating, 
the expansion is faster than in the radiation-dominated era; freeze out occurs earlier, enhancing the X abundance 
at T — Tp [see Eqs. ( |42| ) and ([43|)]. However, as the universe cools from Tp to Trh, the expansion dilutes nx by a 
factor (Trh /Tp) s ; the dilution is more effective than in the matter-dominated case (nx ~ T 3 ) because of entropy 
release during reheating. This explains why Qx is roug hly T RH T F ld I '(T| GW ) 4 times the relic density obtained in the 
case of large reheat temperature. 
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FIG. 1. Shown in the upper graph is the evolution of the X density in the case that the cross section is sufficiently large to 
establish chemical equilibrium prior to freeze out. The lower graph illustrates the case where the cross section is too small to 
establish chemical equilibrium. The two cross sections were chosen to result in the same final X abundances necessary to give 
a critical density of X particles today. In the calculations g f was kept constant at g f = 30. 



An illustration of the freeze out of the X-abundance in equilibrium and out of equilibrium is illustrated in Fig. 



In both cases the final X-abundance is the same. In the top graph the cross section is large enough to establish 
equilibrium prior to freeze out, while in lower graph the cross section is too small to establish equilibrium. 



C. Nonrelativistic Production and Freeze Out 



The relic density calculations performed in Sees. Ill A and IIIB, under the assumptions of out-of-equilibrium and 



equilibrium respectively, match well in the intermediate region in spite of the fact that they are derived with different 
approximations. Indeed, if we use in Eq. ( [45|) the minimum allowed cross section (a s — a s ), corresponding to the 
maximum allowed freeze-out temperature (xf = 5/2), we obtain the bound 
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When this bound is saturated, we are approaching the transition from the results of Sec. Ill E| to those of Sec. Ill A. 
Indeed, the results in Eqs. ( j39|) and ([471 ) turn out to be in fair agreement with each other. Similar conclusions can 
be obtained in the case of p-wave annihilation, but for simplicity in this section we will consider only the case of 
dominant s-wave annihilation. 
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FIG. 2. The shaded areas show the cosmologically excluded regions for a particle of mass Mx with 2 degrees of freedom 
which annihilates in s-wave with a thermal-averaged nonrelativistic cross section {av). The upper-left figure is the usual case 
where particle freeze out occurs when the universe is radiation dominated. In the other frames, we have chosen Mx/Trh = 50, 
100, and 200. The interesting region for cold dark matter (0.025 < Qxh 2 < 1) is between the dashed line and the shaded area. 
The upper right-hand corner of the Mx — (cv) plane is excluded by unitarity arguments. 



In the ordinary case of large reheat temperature {Trh ~ Mx, 



production and freeze out in a radiation- 



dominated universe), Clxh is proportional to 



as can be seen from Eq. (Hq). So except for a logarithmic 



correction, the is no explicit mass dependence to Vlxh- . The constraint from the age of the universe, Vlxh ~ 1, 



implies a lower bound on the X annihilation cross section, as shown in the upper-left graph in Fig. ^|. The unitarity 
limit to (av) as a function of the mass, (ctv)max = 8tt/M x , is also shown. The two bounds cross at the value 
Mx = 340 TeV [|ll], which is usually taken as a cosmological upper bound on any stable massive particle. 

In Fig. U we show how the cosmological bounds on Mx can be relaxed for sufficiently low Trh ■ As Trh is 
decreased, the allowed region in the parameter space (av) versus Mx grows. The numerical results presented in the 
figure (assuming — 30) are in excellent agreement with the analytic estimates. The results are simple to understand. 



In the case of low reheat temperature, the constraint flxh < 1 does not simply give a lower bound on 



but 



rules out a range of val ues. Indeed, for very small (av) we are in the limit of nonrclativistic nonequilibrium production 
discussed in Sec. [II A , and Eq. (|3g) applies. For fixed Trh/Mx, flxh 2 oc M x (av), and the value of (av) forming 



contours of constant Q,xh 2 will scale as M x . As (av) grows for fixed Mx, fix increases and eventually may conflict 
with the age of the universe constraint. After the maximum value of fix is re ached , see Eq. (|39|), a further increase 
of (av) brings us to the limit of nonrelativistic equilibrium production of Sec. Ill B . Now f2x is obtained from Eq. 
(|45|), and the value of (av) to give contours of fixed flxh 2 increases (albeit slowly) as Mx grows. 

Figure || also shows the parameter region in which the particle X could be an interesting cold dark matter candidate, 
0.025 < ilxh 2 < 1. Values of annihilation cross sections and masses which are ordinarily excluded in the case of large 
Trh can now be of particular cosmological and observational interest. 

In Fig. |H we show the cosmologically excluded region (flxh 2 > 1) and the region relevant for dark matter (0.025 < 
Q,xh 2 < 1) as a function of Trh for a fixed value Mx = 100 GeV. Notice how the lower limit on (av), which is 
2 x 10~ 10 GeV -2 in the ordinary radiation-dominated cosmology, is relaxed as Trh is lowered. 
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FIG. 3. The shaded region shows the cosmologically excluded region, as a function of the reheat temperature Trh, for a 
particle of mass Mx = 100 GeV with 2 degrees of freedom which annihilates with a thermal-averaged nonrelativistic s-wave 
cross section (av). The region interesting for cold dark matter (0.025 < Q,xh 2 < 1) is delimited between the dashed line and 
the shaded area. 



A striking implication of dark matter production during reheating is that the unitarity bound, Mx < 340 TeV, 
disappears and one could conceive thermal relics of very heavy particles without conflicting with the age of the 



universe. The necessary assumption is that Tmax is larger than the temperature at which the relevant physical 
processes occur, i.e., T* in the out-of-cquilibrium case and Tp in the equilibrium case. Once this assumption is made, 
the final result on fix does not depend Tmax or other any initial conditions of the inflationary model, but only on 
Trh- Figure || shows how the unitarity bound is modified in presence of low Trh ■ Here we have taken a s = Sir and 
plotted contours of various values of Vlxh 2 as a function of Mx and Trh ■ 




D. Relativistic X 



In this section we consider the case in which the relevant physical processes of particle production and freeze out 
occur when X is still relativistic, and therefore we assume that the inequality of Eq. (|2l|) is not satisfied. If the 
annihilation cross section is not suppressed by any mass scale larger than Mx, then (cry) ~ T~ 2 and X remains in 
thermal equilibrium until it becomes nonrelativistic. More interesting is the case in which the annihilation process 
depends on a new mass scale Mq ^> Mx- Therefore, we define 

rpn 

M - ^ (48) 

for a generic exponent n. As concrete examples, one can think of a heavy neutrino or a neutralino in the regime in 
which the temperature is larger than their masses: in this case n — 2 and Mq is roughly the mass of an intermediate 
gauge boson or of a slepton, respectively. Another interesting example is the gravitino, for which n = and Mq ~ Mpi . 
Finally, one can consider the graviton Kaluza-Klein excitations in theories with 5 large extra dimensions , for which 
n = S and Mq is of the order of the fundamental gravitational scale. 

Let us first consider the case in which X does not reach an equilibrium density. The analysis is similar to the one 



performed in Sec. HI A. At early times, the Boltzmann equation can be approximated by 
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where = g£(3) for bosons, = (3/4)g£(3) for fermions, and g is the number of X degrees of freedom. Using the 
relation between scale factor and temperature in Eq. (jig), we can rewrite Eq. ([l9|) as 
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For n < 6 the A"-particle production dominantly occurs at the lowest possible temperature. In all interesting situations 
we know of, the annihilation cross section is such that n < 6, and therefore we consider only this case. Integrating 
Eq. © up to a final temperature Tf, we obtain for n < 6, 
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If Mx > Trh, Eq. (51) should be evaluated at Tf — Mx and the result be used as an initial condition for the 
nonrelativistic analysis. Since we have linearized the differential equation, within our approximation this simply 
amounts to adding the relic density obtained from Eq. ( |5~T| ) to the contribution derived in Sec. Ill A. 

If Mx < Trh, as it is usually the case for the gravitino, then in Eq. ( [5l| ) we can take Tf — Trh- At temperatures 
smaller than Trh, the universe is radiation dominated and a calculation analogous to the one that led us to Eq. 
( pl| ) shows that in this case A"-particle production dominantly occurs at the largest possible temperature, as long as 
n > — 1. Therefore all the relevant dynamics occurs at T = Trh. The relic a bunda nce can be obtained by rescaling 
the X number density in Eq. (|5l|) to the present temperature, as done in Sec. [II A, to yield 
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The result in Eq. ( J52| ) is valid as long as the X particles do not reach equilibrium at temperatures larger than T RR , 
or n x (T RH ) < n e ^(T RH ). This implies 
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This condition holds for both gravitinos and Kaluza-Klein gravitons. For a Majorana fermion with n = 2, it requires 
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If the condition (|53| ) is not satisfied, then the X particle density thermalizes. However, in this case, they do not 
freeze out before reheating. Indeed, let us consider the ratio of the interaction rate versus the expansion rate at the 
temperature T 
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(56) 



Since for n < 6 the right-hand side is of order unity, Eq. (|56|) shows that if a relativistic X reaches thermal equilibrium, 
then it does not freeze out before reheating. 

In conclusion, for relativistic particles the relevant processes determining their relic abundances occur at reheating 
or afterwards. This, in particular, is true for gravitinos and for the graviton Kaluza-Klein excitations with S < 6 for 
which the cosmological bound derived in Ref. H applies. 



E. Summary of the Different Cases 



Because of all the various cases encountered, it is probably useful to summarize the different possibilities. Let us 
consider a stable weakly-interacting particle X with mass Mx and dominant s-wave annihilation. 

When Mx ~ 17Tmax/4 we are in the deep nonrelativistic regime, and the X relic abundance is strongly suppressed 
by an exponential factor. 

For 17Tmax/4 ^ Mx ~ Tmax we are in the nonrelativistic case. Depending on the value of the annihilation cross 
section, X may or may not reach an equilibrium distribution before freezing out. In the first case (a s > a s ), fix is 
given by Eq. (45), and in the second one (a s < a s ), it is given by Eq. ([56]). The agreement between analytic estimates 
and numerical integration of the Boltzmann equations is illustrated in Fig. |5|. 




FIG. 5. Comparison of numerical vs. analytic results. The nonequilibrium calculation is relevant for a < a, shown in the 
figure. The equilibrium calculation assumes decoupling while nonrelativistic, or xp > 1. 



For Tmax ~ Mx ~ XfTrh, where xf is given in Eq. the particle X is first relativistic, then becomes 

nonrelativistic and finally decouples before reheating. The relic abundance is again given by Eq. (45). 

Lighter X particles thermalize after reheating, erasing any previous information on their number density. The relic 
abundance has the ordinary expr ession given in Eq. (|46[). 

We have also discussed in Sec. HID the case of very light particles with annihilation cross section suppressed by a 
heavy mass scale. In all cases of interest, the relic abundance is determined by the physics at Trh- 



IV. APPLICATIONS TO DARK MATTER CANDIDATES 



A. Supersymmetry 

The neutralino is the most natural cold dark matter candidate in the context of supersymmctric extensions of the 
Standard Model. If the neutralino is dominantly a Higgsino its relic density is typically small, because of the efficient 
coannihilation with other neutralinos and charginos which turn out to be almost degenerate in mass. Upcoming 
LEP2 runs will be able to probe the small window left unexplored in which a light Higgsino could give a significant 
contribution to the present energy density of the universe [£3|. Moreover, if the Higgsino is heavier than the gauge 
bosons, the annihilation channels into and Z° strongly deplete its relic abundance. Significant contributions to 
O then require a lightest supersymmetric particle (LSP) heavier than about 500 GeV, weakening the motivations for 
low-energy supersymmetry. 

The case of a mainly B-ino lightest neutralino is much more promising for dark matter. First of all, we should recall 
that most of the supersymmetric models obtained from supergravity usually predict that the Higgs mixing parameter 
fx is large. This is because /i should compensate the large radiative corrections to the Higgs mass parameters in order 



to achieve the correct size of electroweak symmetry breaking. Therefore, a common expectation is that the lightest 
neutralino is an almost pure i?-ino. 

In the early universe, the £?-ino will mainly annihilate into fermion pairs through ^-channel exchange of squarks 
and sleptons. Exceptions occur only for pathological situations in which there is a resonant s-channel exchange of Z° 
or a Higgs boson. Actually, because of the large hypercharge of the right-handed electron and the expected lightness 
of sleptons compared to squarks, it is often a good approximation to include in the annihilation cross section only 
the exchange of the right-handed sleptons. Summing over three slepton degenerate families with mass rhe R , the £>-ino 
annihilation cross section parameters are 

247ra 2 / \ ~ , s 

«°=°> ^=c^( l+ 4J • (57) 

Notice that the annihilation process is p-wave suppressed because of the Majorana nature of the neutralino. 

Cosmological considerations give an upper bound to the B-ino mass. Indeed, the requirement that charged particles 
are not the LSP implies rhe R > Mb- The minimum allowed B-ino relic abundance corresponds to the maximum 
annihilation cross section and therefore to the minimum rht R . Setting fht R = Mb in the expression for f2, one obtains 
an upper bound on the i?-ino mass of about 300 GeV (for £lh 2 < 0.3) |05). The bound can be weakened in the presence 
of resonant s-channcl annihilations, once a small Higgsino admixture is introduced. Moreover, in constrained models 
in which the supersymmetry-breaking masses satisfy simple universal relations at the GUT scale, this bound reduces 
to 200 GeV Q. 

However, as emphasized in Ref. |jT5| , whenever the sleptons and the .B-ino become degenerate in mass within 
about 10-20%, one cannot ignore the effects of coannihilation. These effects can modify significantly the B-mo relic 
abundance, because annihilation channels involving the charged sleptons have large cross sections which are not p- 
wave suppressed. Indeed, even in the case of the constrained model, the previous limit on the S-ino mass can be 
relaxed to about 600 GeV 
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FIG. 6. The allowed region in the MB-rhi R plane in the standard cosmology where freeze out occurs in the radia- 
tion-dominated era is indicated by the shaded region above the disallowed region where the LSP is charged and below the 
curve marked "forbidden Mb/Trh — > 0" above which S)g/i 2 > 1. Coannihilation effects, here neglected, modify the bounds in 
a narrow region where B and £r are nearly mass-degenerate. If the B-ino freezes out during reheating the forbidden region 
where fish 2 > 1 is a banana-shaped region. The allowed region of parameter space is above the charged-dark-matter region 
and outside the banana-shaped region. The size and location of the banana-shaped disallowed region depends on Mb/Trh- 
Some examples are shown in the figure. 



Coannihilation effects do not significantly modify the bound on the slepton mass for a fixed value of Mb (as long 
as it is not too close to fht R ). On the other hand, these bounds can rapidly disappear if we consider low values 
of the reheat temperature. The effect is illustrated in Fig. 0, which shows the values of slepton and B-ino masses 



incompatible with the constraint flh 2 < 1, for different choices of Trh, in the case of a 100% pure -B-ino LSP. The 
upper bounds on fhi R are drastically relaxed.^ Therefore, this shows that cosmological arguments based on relic 
abundances, used to set upper bounds on supersymmetric particles rely on specific untested assumptions. A low 
reheat temperature can completely change the picture. 

The same arguments can be applied to other supersymmetric dark matter candidates. One possibility which one 
encounters in theories with gauge-mediated supersymmetry breaking |25| , ^6|| is given by the messenger scalar particle 
with the same gauge quantum numbers of the neutrino. The cosmological upper bound on this particle mass is about 
3 TeV, much lower than the natural theoretical expectation. A low reheat temperature easily relaxes the bound. 

The same can be said about the limits on unstable particles whose decays into energetic products may jeopardize 
the successful predictions of nucleosynthesis Jig] . It is clear that our results relax the bounds in the {Mx,tx) plane, 
where Tx is the decay lifetime of the unstable particle. 



B. Massive Neutrinos 



Let us now analyze the implications of our approach for massive neutrinos, in the case in which they are stable. 
First, let us briefly recall the standard (i.e., when the reheat temperature is very large) prediction for the abundance of 
massive neutrinos. Neutrinos are initially kept in equilibrium by weak interactions. For neutrinos lighter than about 
1 MeV, freeze out occurs at To — 2.3 MeV for electron neutrinos and To ~ 3.7 MeV for muon and tau neutrinos, 
so the neutrinos are relativistic at freeze out. The current abundance of a generic relativistic particle can be easily 
estimated to be || 

2 g x 10.75 m x fKQ s 
Qxk ^-yjT^MeV- (58) 

Here gx is the number of degrees of freedom of the species. Standard Model (SM) relativistic neutrinos decouple from 
chemical equilibrium when <?*(Td) ~ 10.75. One usually concludes that the mass of SM neutrinos cannot be larger 
than 94/i 2 eV, the Cowsik-McClelland bound @. 

If neutrino masses are larger than the freeze-out temperature, they decouple from the thermal bath when they are 
nonrelativistic. In this case the annihilation cross section is proportional to G^to 2 , and requiring Q, v h? < 1 provides 
an lower bound on m v of about 2 GeV, the Lee- Weinberg bound |l7j . This means that neutrino masses in the range 
94 eV ^ m v ^ 2 GeV are cosmologically ruled out. This is the celebrated Cowsik-McClelland-Lee- Weinberg (CMLW) 
bound ^60. 

This picture has to be modified in the case in which the reheat temperature is small. The standard CMLW bound 
is based on the assumption that neutrinos have reached thermal and chemical equilibrium in the radiation-dominated 
universe. This is equivalent to the hypothesis that the maximum temperature obtained during the (last) radiation- 
dominated era, that is, the reheat temperature Trh , is much larger than the decoupling temperature To- We have 
no physical evidence of the radiation-dominated era before the epoch of nucleosynthesis {i.e., temperatures above 
about 1 MeV). Therefore, let us explore the possibility that the largest temperature of the Universe during the 
radiation-dominated phase is very small. Indeed, it has been recently shown that the smallest value not excluded by 
nucleosynthesis at more than 95% CL is T RH ~ 0.7 MeV 

Since neutrinos have only weak interactions, it is very difficult for the thermal scatterings of particles during 
the reheat stage to generate SM neutrinos through processes like e + e~ — > vv and to bring neutrinos into chemical 
equilibrium. Furthermore, decreasing the reheat temperature increases the rate of the expansion of the Universe, as 
explicitly seen in Eq. (pp|), making it more difficult for the weak interactions to bring the neutrinos into chemical 
equilibrium. Therefore, if the reheat temperature is small enough, one should expect that the SM neutrinos produced 
during the reheat stage never reach chemical equilibrium. In other words, at the beginning of the radiation-dominated 
phase neutrinos populate the thermal bath, but they have a number density, n Vl which is smaller than the equilibrium 
number density, n e ^ . This is well- illustrated in Fig. |t], which shows the evolution of the number density of a muon 
or tau Dirac neutrino with mass 12 keV, for Trh = 1 MeV. It is clear that the number density only grows, and 



5 This is particularly welcome in those scenarios where the supersymmetric flavor and CP problems are avoided if the first 
two generations of sfermions are heavier than a few TeV and approximately degenerate in mass. If the lightest supersymmetric 
particle is essentially _B-ino-like then requiring that all flavor changing neutral current and CP-violating processes are adequately 
suppressed imposes a lower limit on the B-ino mass of typically 200 to 300 GeV VM. 



it is always too small to catch up with the equilibrium number density, i.e., n v -C n\ 
annihilations arc not efficient. 



This implies that neutrino 
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FIG. 7. The evolution of the number density per comoving volume of Dirac neutrinos for m v 
This figure illustrates the fact that neutrinos never attain chemical equilibrium. 



12 keV and T RH = 1 MeV. 



This result applies both to relativistic and nonrelativistic SM neutrinos, and implies that the present abundance of 
neutrinos in low Trh models is much smaller than predicted assuming that the largest temperature of the radiation- 
dominated Universe was much larger than a few MeV. This is the reason why the CMLW bound on neutrino masses 
is significantly relaxed in low Trh models. 

In Fig. ^ we present a full numerical computation of the abundance of massive neutrinos in terms of the total 0„/i 2 
as a function of the mass of the neutrino, for different values of the reheat temperature. In all our results we have 
numerically solved the Boltzmann equation making use of the exact definition of the thermally averaged cross section 

M 
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AmlTKl{m v /T) J 4m2 



/ ds avE u E B y/a - Ami Ki{V~s/T) , 

J 4m? 



(59) 



where s is the center-of-mass energy, Ki are the modified Bessel functions and we borrowed the expressions for avE^Ep 
from the appendix of Ref. .[] 

For a Dirac neutrino, the treatment of the number of degrees of freedom requires some attention. In the relativistic 
case, only two of the four helicity degrees of freedom are produced because the generation of the wrong-helicity states 
is suppressed by (m^/2E^) 2 . On the other hand, in the nonrelativistic regime, all four degrees of freedom interact 
with full strength and will be produced during the reheating stage. Here, we make the "helicity approximation" and 
assume two degrees of freedom for a relativistic neutrino species and four for a nonrelativistic neutrino species. 

At low reheat temperatures neutrinos of a given family v a can be produced by the processes e + e~ — > v a v a and 
vpvp — > v a v a , where (3 is different from a. The inclusion of the e + e _ scattering in the integrated Boltzmann equation 
for the number density of neutrinos n Va is straightforward because e~ and e + are kept in chemical equilibrium by 
the fast electromagnetic processes and Eq. (||) applies. However, the inclusion of the vpv$ scatterings is more delicate 
because vp and z/g are themselves not in chemical equilibrium. A complete solution of the problem would require a 
detailed kinetic treatment of all the neutrino distribution functions f v (p,t) in momentum space. This computation 
is now in progress. In this paper we have limited ourselves to include the vpvp scatterings by defining in Eq. (59) an 



effective cross section a 
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. To convince oneself of the validity of this approximation, 



one may notice that we recover the usual standard CMLW bound when the reheat temperature is larger than about 
7 MeV. This is in agreement with the results obtained in Ref. p7[| where the Boltzmann equations in momentum 



6 In Ref. [B9|, for the Majorana case, CV and Ca have to be interchanged. We thank G. Raffelt for communications about 
this point. 



space were numerically solved for massless neutrinos and it was shown that for Trjj ;$ 7 MeV the effective number of 
massless neutrinos N v = p v j p^ q starts deviating from 3. 
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FIG. 8. The contribution to the closure density of Dirac and Majorana neutrinos for different values of the reheat temperature 
(g, — 10.75 is used in the definition of Trr)- 

From Fig. ^|we infer that in the case in which neutrinos are relativistic, Q u h 2 is approximately given by 

3 



£l Utl h 2 = VL VT h 2 ~ ( 12 Y y ) ( mv ) relativistic Majorana and Dirac 



(60) 



The full numerical calculation shows that neutrino masses as large as m v ~ 33 keV are compatible with Q, v h 2 < 1 
for the limiting reheat temperature Trh — 0.7 MeV. This shows that SM neutrinos with masses up to about 33 keV 
are perfectly compatible with cosmology and may even play the role of warm dark matter. 

Let us now briefly see what our findings are when the neutrinos are heavier than an MeV. This is possible only 
for the tau neutrino, for which the present experimental upper limit on the mass is 18.2 MeV [ j30| . For m VT ~ 10 
MeV the ordinary calculation of the relic abundance in a radiation-dominated universe predicts a value of fi„ T ~ 10 4 , 
which is definitely excluded. Our findings indicate that 



o ^ ( T Rh\ 7 /14MeV\ 3 

ll v h = nonreiativistic Dirac , 



n Ur h 2 = 



(TrhV ( 13 MeV 
I MeV J 



nonreiativistic Majorana . 



(61) 



The strong dependence upon th e rehe at temperature is easily understood by realizing that for Trh hi the MeV range 
we are in the condition of Sec. IIIB| (a s < a s and a p < a p ), and therefore the neutrino relic abundance is given 



approximately by Eq. ([36]). Our numerical result is very well reproduced analytically by taking Eqs. ([27]) and ( |3(: 
and using, for a Dirac neutrino, 



G F m Vr ( 1 _ 2 sin2 &w + 4 gin 4 Qw j ^ ( 62 ^ 



and, for a Majorana neutrino, 



2?r \2 



2(,/ '" v ^i-2sin 2 6» W /+4sin 4 6»vi/ I . (63) 



These values are obtained taking into account only the process e + e — > v r v r in cr e // because, as we have numerically 
checked, the contribution to the production of heavy tau-neutrinos from light neutrino annihilations is negligible at 



the production time, see Eqs. (29) and (pOp. 

It is interesting to notice that there exists a small window of m„ T and Trh for which the ordinary tau neutrino is 
an a ccep table candidate for cold dark matter. However, because of the large powers of m VT and especially of Trh in 
Eq. (|6l|), the allowed window is very limited. Of course in this scenario there must be a conserved quantum number 
to keep the tau neutrino stable in order to be the dark matter; this hypothesis is in apparent contradiction with the 
experimental evidence on neutrino oscillations. 

It has also been proposed that the tau neutrino could be a cold dark matter candidate if it had a magnetic moment 
of the order of 10 -6 Bohr magnetons pH]] . The present experimental limit on the v T magnetic moment of 5.4 x 10~ 7 /is 
rules out this possibility. We find that choosing Trh as low as 1 MeV can rescue this possibility, since f2„ T /i 2 ~ 1 for 
a tau neutrino magnetic moment with the maximum allowed value. 

We conclude that the cosmologically disallowed region for neutrino masses is 

33 keV < m v < 6 MeV (64) 

for a Dirac neutrino and 

33 keV < m v < 5 MeV (65) 

for a Majorana neutrino. We stress that so long as the reheat temperature is an unknown parameter, this should 
be considered the real CMLW bound on neutrino masses. Our findings indicate that neutrinos can still play the 
role of warm or cold dark matter and that the impact of massive neutrinos on nucleosynthesis has to be revisited. 
These and other issues are currently under investigation. Finally, we note that the above values were found assuming 
Maxwell-Boltzmann statistics and assuming the annihilation products are in equilibrium. 



C. Cosmological Bound on Axions 

The invisible axion is still the most elegant solution to the strong CP problem |32|,[33). An axion model has one 
basic free parameter, the axion mass m a , or equivalently, the Peccei-Quinn (PQ) symmetry breaking scale Jpq. The 
mass and symmetry breaking scale are related by 

mQ " - 62eV l^' (66) 

where N is the color anomaly of the PQ symmetry. 

Several astrophysical lower limits on fpQ are based on the requirement that the axionic energy losses from stars, 
notably red-giant stars, globular-cluster stars, or the core of supernova 1987A, are not in conflict with the observed 
properties of these objects (34). These limits imply m a ^ 10 -2 eV (equivalently Jpq > 10 9 GeV), indicating that 
axions, if they exist, are both extremely light and very weakly interacting. 

An upper bound on the PQ scale comes from cosmological considerations. Let us call O the strong CP-violating 
phase. Today O is anchored at the CP-conserving value, O = 0. However, the axion mass is very temperature- 
dependent |35[ 

m a (T)c 0.1 m a (^p) 3 7 , (67) 



where Aqcd — 200 MeV is the QCD scale and the relation is valid for T ^> Aqod/t- At very high temperatures 
the axion is essentially massless. This means that no special value of is specified by the dynamics, and all values 
of the phase are equivalent. The axion mass turns on at a temperature T\ such that m a (Ti) ~ 3H(Ti), and the 
axion field starts evolving toward the minimum at = 0, eventually oscillating around it. These cosmic oscillations 
of the axion field correspond to a zero-momentum condensate of axions which does not decay. The energy density in 
axions today from this misalignment mechanism exceeds the critical density unless fpQ < 10 12 GeV pq| . This result 
has always been considered particularly disappointing from the theoretical point of view, since weakly coupled string 
theory possesses numerous axion candidates, whose decay constant is, however, of the order of the string scale and 
therefore much larger than 10 12 GeV. 

The purpose of the present section is to demonstrate that the cosmological bound on the PQ scale is significantly 
relaxed if we make the assumption that the reheat temperature is smaller than the QCD scale (similar considerations 
have also been made in Refs. ]37| , |38| ] ) . 

Suppose then that Trh < Aqcd ~ Tmax- This means that the axion coherent oscillations commence when the 
universe is still matter-dominated and reheating is not completed. During this epoch the Hubble rate is given by Eq. 
©, and the axion mass is still given by Eq. ( p7| ) since the universe is populated by a thermal bath between Tmax 
and Trh- Axions start oscillating at a temperature T\ when m a (Ti) ~ 3-ff(Ti): 
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(68) 



At temperatures T < T\ the number density of axions scales like a 3 even though the axion mass is still varying. At 
the reheat temperature we have 



n a (T RH ) = n a (Ti) 



a{T RH ) 



( fp Q \ 2 9 3 * /2 (Trh) T% h 
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<7„(Ti) TfM P i 



(69) 



where 0i is the initial displacement of the CP-phase. Notice, in particular, that the ratio n a /s does not remain 
constant during the cosmological evolution from the temperature T\ to Trh- This is because there is a continuous 
release of entropy. However, when reheating is completed and the universe enters a radiation-dominated phase, the 
ratio n a /s is conserved, and one can easily compute the present abundance of axions from the misalignment mechanism 
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Requiring that Q a h 2 < 1 gives 
fpQ 
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(70) 



(71) 



Therefore, the cosmological axion problem is ameliorated]^] Furthermore, in the strong coupling vacuum described in 
Ref. |39|, the QCD axion might be a boundary modulus. Dimensional analysis suggests approximately 10 16 GeV for 
the decay constant of such a boundary axion, not necessarily in contradiction with the upper bound of Eq. ( |7l| ) , see 
also the discussion in Ref. |f40| . 

V. IMPLICATIONS FOR GUT BARYOGENESIS AND LEPTOGENESIS 

The explanation of the observed baryon asymmetry (B) in the early universe, of the order of 10~ n in units of 
the entropy density, remains a fundamental cosmological question pl|. Several theories with typical energy scale 



7 In cosmologies where the universe is dominated early on by the coherent oscillations of some moduli field the axion bound 
is significantly weakened, as had been already observed in Ref. J^Q. In Ref. j38), however, the estimate on the upper bound 
on fpQ did not take into account the presence of the thermal bath before the completion of reheating, and therefore neglected 
the dependence of the axion mass on the temperature. 



much higher than the electroweak scale can explain the observed baryon asymmetry. For instance, in Grand Unified 
Theories (GUTs) the out-of-equilibrium decay of heavy Higgs particles may be responsible for the direct generation 
of the baryon asymmetry |42] ]. Alternatively, the baryon asymmetry may be produced from a lepton asymmetry 
(L) ^] using the fact that any lepton asymmetry is reprocessed into baryon number by the anomalous sphaleron 
transitions p4[ . In the simplest scenario, the lepton asymmetry is generated by the out-of-equilibrium decay of a 
massive right-handed Majorana neutrino, whose addition to the Standard Model spectrum breaks B — L. 

However, any scenario for the generation of the baryon asymmetry based on the out-of equilibrium decay of some 
heavy particle X depends crucially on the assumption that these particles were nearly as abundant as photons at 
very high temperatures. This imposes a lower bound on the reheat temperature, Trh ^ Mx- On the other hand, in 
supersymmetric models the requirement that not too many gravitinos are thermally produced after inflation provides 
a stringent upper bound on the reheat temperature of about 10 8 to 10 10 GeV ||. If this bound is violated, the decay 
products of the gravitino destroy light nuclei by photodissociation and hadronic showers, thus ruining the successful 
predictions of nucleosynthesis. Therefore, any out-of-equilibrium decay scenario would require Mx ~ 10 8 -10 10 GeV, 
a condition which looks particularly problematic for GUT-inspired baryogenesis. 

In order to relax this limit one usually envisages two possibilities. Either, the heavy particles are produced directly 
through the inflaton perturbative decay process |^5| (this requires that the mass of the inflaton is larger than Mx) 
or they are generated through nonperturbative process taking place at the preheating stage (see Ref. (Hw in the case 
of GUT Higgs boson induced baryogenesis and Ref. |47| in the case of leptogenesis). 

In this section we wish to show that the heavy particles X may be abundantly produced by thermal scatterings 
during the reheat stage even though the reheat temperature Trh is smaller than Mx ■ Again, this is made possible by 
the fact that Trh is not the maximum temperature of the universe during reheating. For the sake of concreteness we 
will focus on the leptogenesis scenario, but our findings can be easily generalized to any out-of-equilibrium scenario. 

Let us indicate by N = n^a 3 the number density per comoving volume of the lightest right-handed neutrino, 
the one whose final decay (into left-handed leptons and Higgs bosons) is responsible for the generation of the lepton 



asymmetry. Following the notations of Sec. II A, we can write the Boltzmann equation for N as 

(72) 



dN _ c N A 1 / 2 (N - N eq ) 



dA 



where 



c N = J^-^ (T Nl + 2 r M + 4 r M ) . (73) 

V 8lT 1 RH 

Here is the decay rate of Ni (for the processes Ni — > H^£l,H£l); ^h,a and Th,t are the rates of the scattering 

— (3) ~ 

processes containing Ni in the final state, mediated by the Higgs boson in the s channel {tR,q L — > IlN\) and in the 

(3) 

I channel (Il<I l —> t^Ni), respectively. 

Let us suppose first that Tmax ~ M\, where M\ is the N\ mass. Under this assumption we have 



a? m r - 3A?A? T3 r - 3A ? A ? Tln M i (74) 

A| --M 1; T Ks -^— r r feii --^-3-Tln— , (74) 



where = (AA T )n, with Xij Yukawa coupling of Ni, and \ t is the top Yukawa coupling. For a more transparent 
physical interpretation, it is convenient to express A^ in terms of the parameter 

mi = 2V2G F M 1 ' 

In the limit of small mixing angles the parameter mi coincides with the mass of one of the light (mainly left-handed) 
neutrinos. 

If the right-handed neutrinos do not reach an equilibrium density (N <C N eq ), we can approximate Eq. ( [72] ) by 

dN cnA 1 ' 2 N eq 

(76) 



dA JW, 



Along the same lines of Sec. Ill A, we can integrate Eq. ( J76| ) by approximating it to a Gaussian integral in the full 
range between ^4=1 and A — oo. This is a good approximation because the exponential suppression in N eq makes 
the right-hand side of Eq. (76) negligible anywhere outside a small interval of scale factors centered around A = A* 



corresponding to T* = Mi/ 10 for the inverse decay, and to T* = Mi/9 for the Higgs-mediated AL = 1 processes. It is 
easy to show that the main source of right-handed neutrinos is represented by the inverse decays whose contribution 
to is given by 

^ c r Nl M 3 Pl HjT* H 

°° -3/2 M? U0 

where 
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c = 



= 2.7 x 10 3 . (78) 



Notice that the final abundance is suppressed only by powers of the right-handed mass, there is no Boltzmann 
suppression exp(— Mx/Trh)- Furthermore, the abundance is proportional to the rate of production (accumulation) 
Tn ± - However, the consistency condition Nqo ^ N eq (T*) gives an upper bound on Tn ± 

r < 77 V2 M l / 70 A 

rjVl ^ 8000 ff * M^Ih' { ) 

or, equivalently, an upper bound on m\ 

< 8 »> 

These bounds can also be expressed in terms of a more familiar quantity 



H 



< 1 . (81) 

T=T» 



This condition assures that when the right-handed neutrinos are produced, their direct decay is inefficient. The 
limiting case ~ 1 would mean that the right-handed neutrinos enter into chemical equilibrium as soon as they are 
generated. 

The right-handed neutrinos may decay before or after the universe reaches the reheat temperature Trh, depending 
on the value of Ai. Suppose that they decay after the end of the reheat stage (which is, for instance, always the case 
when M\ < 8Trh)- This means that at T = Trm, the ratio of the number density of right-handed neutrinos to the 
entropy density is given by 



n Nl _ 257TC T Nl MpiT 7 RH 



s{Tr H ) M? 

This ratio remains constant until the right-handed neutrinos decay generating a lepton asymmetry 



(82) 



10 5 e Y Nl M Pl Tl H 

~ JT2 Wl ' (83) 

y* 1 

where we have indicated by e the small parameter containing the information about the CP-violating phases and 
the loop factors and we have again taken into account the factor of 1/8 due to the release of entropy after Trh- 
The corresponding baryon asymmetry is B = (28/79)L, assuming only Standard Model degrees of freedom Q], and 
therefore 

By virtue of the bound of Eq. (f79|), this baryon asymmetry is constrained to be smaller than the critical value of 
2(e/<?*)(Tfl/f/Mi) 5 (and of 2 x 10~ 5 e/<7* if we use the constraint T» > Trh)- The requirement that B is larger than 
2 x 10~ n implies 



It is easy to convince oneself that this is also the result in the case in which the right-handed neutrinos decay before 
the reheat stage is over. Equation (|8^) provides a necessary condition on the mass Mi of the lightest right-handed 
neutrino in a leptogenesis scenario, correcting the naive estimate Mi < T R h- The relaxation of the naive bound by 
more than one order of magnitude is certainly welcome to make leptogenesis more compatible with the cosmological 
gravitino problem. 

Let us suppose now that Tmax ~ Mi > Trjj, and that inverse decays or production processes containing the Ni in 
the final states can bring the right-handed neutrinos to equilibrium before they become nonrelativistic. This amounts 
to requiring that the AL = 1 interactions with total rate 7^ = (T^x + 2 Th, s + 4 T^j) are in thermal equilibrium at 
T ^ M\. Therefore the standard out-of-equilibrium parameter 

K^ 1 -^- (86) 

ti T=Mi 

is larger than unity. 

In this case, the lepton asymmetry can be written as 

45 ez 3 / 2 e~ z t fT RH \ 5 

L = 2 9/2^/2 - {-Tf) ■ ( 87 ) 

Here Zf = Mi/Tf, where Tf is the temperature at which the processes that damp the baryon asymmetry go out of 
equilibrium, and the last factor in Eq. (|87| ) accounts for the dilution caused by the expansion in the pre-reheat phase. 
We are assuming that Tf is larger than the reheat temperature Tret. 

If the inverse decay dominates over scatterings (we will later quantify this condition), then Tf is approximately 
determined by 

Tin = H\ T=T , (88) 

where the inverse decay rate at T < Mx is 

Fid = ^zV 2 e-*T Nl . (89) 



Here z = Mi/T and is given in Eq. ([74]). In terms of the parameter K, defined by 
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(90) 



becomes KzV~' 2 e Zf ~ 1, which is approximately solved (for K < 10 5 ) by Zf ~ 16X 06 . Replacing Zf in 



Eq. 

Eq. 087|), we find 



45 (T RH \° e 
V2^/2 ^ Mi J ff .JS-0-M ■ iy J 



Finally, the baryon asymmetry B — (28/79)L is given by 



100\"Vl(r 3 eV\ fT R „\\ : ^ n _ 3 



Equation (|9^) is valid as long as inverse decay processes dominate over AL = 2 scattering processes in damping 
the baryon asymmetry. Let us now study the condition under which this hypothesis is justified. The rate for lcpton- 
violating scatterings mediated by Ni exchange in the sort channels, for T < Mi, is given by 



AL=2 



lG 2 F Mlm\i 
2tt 3 



(93) 



Therefore, the AL = 2 scatterings are out-of-equilibrium at high temperatures and equilibrate at a temperature 
corresponding to 



zF=* ^^Mi 



The assumption that led to Eq. ( |92| ) is then valid as long as z^ L ~ 2 > Mi /Trh, which implies 

/ g, \ 1/4 A0 10 GeV\ 1/2 ^ 

TOl< (loo) N^J cV - (95) 

When the condition (|9^) is not satisfied, AL = 2 scatterings lead to an exponential suppression of the baryon 
asymmetry. 

The baryon asymmetry in Eq. (|92| ) is valid when three conditions are verified: the right-handed neutrinos N% 
reach equilibrium, which implies K > 1; the temperature Tf is larger than Trh, which implies Mi > f6 Trh', 
Eq. ( p5| ) is satisfied. These three conditions together imply a maximum value of Mi = (I00/ 5 *) 1/2 3 x 10 11 GeV, 
corresponding to a maximum 

T RH = (100/ff*) 1 / 2 2 x 10 10 GeV. The maximum baryon asymmetry, achieved when 
K ~ f , is £? ~ e(f00/,g*)(Tfjjy/Mi) 5 2 x I0~ 3 . Therefore, in this case, a sufficiently large baryon asymmetry requires 
Mi ~ 10 T fl jj and e ~ f 0~ 3 . 



VI. CONCLUSIONS 



In this paper we studied the observational consequences of having a reheat temperature Trh smaller that the 
characteristic temperature at which a certain cosmological process occurs. We first described the dynamics of reheating 
and derived general expressions for the relic abundance of particles whose standard freeze-out temperature is larger 
than Trh- For nonrelativistic particles we found two different regimes. If the annihilation cross section (av) is smaller 
than the critical value in Eqs. (|37|-|3S|), than the present relic abundance is proportional {av), see Eq. (p6|) . In the 
other case, the relic abundance is inversely proportional to {av), sec Eq. (|45|), but because of the fast expansion before 
reheating its expression differs from the usual result in a radiation-dominated universe. 

We applied our general results on relic abundances in \ow-Trh cosmologies to several cases of interest. A first result 
is that the usual unitarity bound of 340 TeV on stable-particle masses can be relaxed. The new excluded ranges of 
stable-particle masses as functions of Trh are shown in Fig. |]. 

We revisited the parameter regions of supersymmetric models leading to viable cold dark-matter candidates in the 
light of low Trh- In particular, we found that the upper bound on the slepton mass, as a function of the LSP B-ino 
mass, can be significantly relaxed, as quantitatively shown in Fig. || Large regions of parameter space that have been 
considered to be ruled out by cosmological arguments can instead give a relic neutralino density close to the critical 
value. 

Next, we considered how a low reheat temperature can reduce the relic abundance of massive neutrinos, here 
assumed to be stable. The requirement that tt u h 2 < I gives the CMLW bound m v < 94 eV only if Trh is larger 
than about 7 MeV. The bound becomes significantly weaker for lower Trh; for instance, m v < 12 keV for Trh = 1 
MeV. For very massive neutrinos the Lee- Weinberg bound, m v > 2 GeV, is also modified. Again for Trh = I MeV, 
the limit becomes m v > 14 (13) MeV for a Dirac (Majorana) neutrino. This means that there is even the possibility 
that a stable v T with mass consistent with the direct experimental limit (m UT < 18.2 MeV) freezes out when it is 
still nonrelativistic and becomes cold dark matter. Indeed, for Trh — 0.7 MeV, which is the lowest value of Trh not 
excluded by nucleosynthesis, the constraint Sl VT h 2 < 1 excludes only the range of masses 33 keV < m„ T < 6 MeV 
(Dirac) and 33 keV < ui Vt < 5 MeV (Majorana). This result resurrects the possibility of neutrinos as warm dark 
matter. For instance, if Trh = 1 MeV, f2„/i 2 = 0.3 for a neutrino of mass 4 keV. 

The requirement that the energy stored today in the axion oscillations (caused by a misalignment between its 
high-temperature and low-temperature configurations) is not larger than the critical value imposes a bound on the 
Peccei-Quinn symmetry breaking scale, fpQ $ I0 12 GeV. If Trh is less than Aqcd, this bound can be relaxed to 
values close to the GUT scale. 

Finally, we have investigated the impact of low Trh in the explanation of the observed baryon asymmetry by a 
leptogenesis mechanism. In this context, low Trh means Trh < Ml, where Mi is the mass of the lightest of the three 
right-handed neutrinos. Therefore, we discuss here values of Trh much larger than in the previous cases. However, 
the formalism is the same, because what matters is that Trh is less than the relevant physical energy scale. We have 
found new expressions of the baryon asymmetry as functions of Trh, see Eqs. ( jgj|) and (p2|). Moreover, values of Mi 
an order of magnitude larger than Trh can still lead to a sufficient density of right-handed neutrinos to explain the 
baryon asymmetry. This is a welcome and important result, when one tries to make leptogenesis consistent with the 
cosmological gravitino problem. 
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APPENDIX: EVOLUTION AFTER T RH 

Numerical results show that there is about a factor of 8 increase in the comoving entropy after Trh, and that the 
value of $ at T = Trh is &rh = 0.79$/, i.e., only about 21% of the comoving <j> energy density has been extracted 
at reheat time. In this Appendix we will give an explanation of these results. 

Assume that at Trh the radiation energy density is less than the <f> energy density. Then the evolution equation 
for $ is 



/ 2 \ 1/2 



AV2$l/2 _ ( A1) 



dA V 30 

Integrating this equation from A = 1 to A = Arh with the initial condition $ = we obtain 



$ 7 ) 3 V 30 

Evaluating Eq. fllS|) at T — Trh we obtain 



= l-;hr ^i rh (Arh* I). (A2) 
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Using this relation, we can estimate 



<Z>RH /13 



$r V 15 



0.75 . (A4) 



This is in good agreement with the numerical results and is independent of model parameters. 

The second step is an estimate of how much entropy is released after Trh, if &rh = 0.75$/. Let us make the crude 
approximation that the entropy release is instantaneous just after Trh- Then ApR = Ap^, and therefore 

Ap R = 0.75 . (A5) 

A RH 



Using Eq. ([A3; ) we obtain 

Ap R = 0.16 7r 2 g*T RH , (A6) 

which gives 

Ap R E 



PR 



(A7) 



Therefore the analytic estimate in the case of the instantaneous approximation predicts a model-independent factor 
of 5 3 / 4 ~ 3.4 for the release of entropy after Trh- The model dependence comes in if we relax the assumption of 
instantaneous release of entropy, increasing the estimate, since p^ redshifts like A~ 3 . The numerical result shows an 
increase of about a factor of 8 in the comoving entropy after Trh- 
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